AD752239 


A  REVIE(7  OP  VIBRATIONAL  CROSS  SECTIONS 
INVOLVING  H2O,  CO2,  H2,  CO  and  He 


r  APCRl-72-0539 


Edward  R.  Fisher 


Wayne  State  University 
Research  Institute  for  Engineering  Sciences 
Department  of  Chemical.  Engineering 
Detroit,  Michigan  48202 

Contract  No.  F  1S528-72-C-0007 


SPECIAL  REPORT 
August  1972 

Contract  Monitor:  Alva  T.  Stair 
Optical  Physics  Laboratory 


Approved  for  public  release;  distribution  unlimited 


Sponsored  by 

Defense  Advanced  Research  Projects  Agency 
ARPA  Order  No.  1856 
Monito’.  ed  by 

Air  Force  Cambridge  Research  Laboratories 
Air  Force  Systems  Command 
United  States  Air  Force 
Bedford,  Massachusetts  01730 

Reproduced  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

U  S  Deportment  of  Commerce 
SprtngfteM  VA  22151 


ARPA  Order  No.  185 S 
Program  Code  No.  2E50 

Contractor;  Wayne  State  University 
Effective  Date  of  Contract;  12  July  1971 


Contract  No.  F19628-72-C-0007 

Principal  Investigator  and  Phone  No^ 
Dr.  Adolf  R.  Hochs tim/31 3  577-3867 

Project  Scientist  and  Phone  No 
Dr.  Alva  T.  Stair,  Jr./617  86i-49li 

Contract  Expiration  Date;  12/31/72  * 


apply  to  the  National  Teohnical  fnfoiilia?Srie?5?Sf 


$ 

- ■-  . 


Unclassified 

Security  CInssification 


.  DOCUMEMT  CONTROL  DATA  -  R&D 

U  /ro,pn' aniotation  muu  be  entered  when  the  overall  report  is  clussif,ed) 

wayne  btate  University 


Research  Inslitute  for  Engineering  Sciences 

^  Engines  A,g 

S-  NEPORT  TITLE  - * - 


2a.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


2i>.  CROUP 


3.  REPORT  TITLE  - - - - 

H^f  mYnd  SECTIONS  INVOLVING  H^O,  CO^, 


4.  DESCRIPTIVE  NOTES  (Type  of  report  and  inclusive  dates) 

Scientific  Interim 


S.  AUTHOR(s)  (First  none,  middle  initial,  last  rume ) 


Edward  R.  Fisher 


«•  REPORT  DATE  '  ~  - - - 

August  1972 

•a.  CONTRACT  OR  GRANT  n'HT'aKjJa  Order  No.  1B56 


F19628-72-C-0007 

b.  PROJECT.  TASK.  WORK  UNIT  NOS. 

8692  -  n/a  -  n/a 

c.  DOD  ELEMENT 

62.3010 

•L  DOD  SUBELEMENT 

Til  a. 


Tt--  TOTAL  NO.  OF  PACE.S 

ao 


So.  ORIGINATOR’S  REPORT  NUMBEt^S; 

Special  Report 


70.  NO.  OF  REFS 
28 


10.  DISTRIBUTION  STATEMENT  - 

A  -  Approved  for  public  release;  distribution  unlimited. 


AFCRL-72-0539 


H.  SUPPLEMENTARY  NOTES 

This  research  was  supported  by  the 
Defense  Advanced  Research  Projects 
Agency 


T£  SPONSORING  MILITARY  ACTIVITY 

Air  Force  Cambridge  Research 
Laboratories  (OP) 

L.  G.  Hanscom  Field 
Bedford,  Massachusetts  01730 


^7 - - - - - - -  1  - '  UltJU 

.  , Theoretical  vibration  cross  sections  calculated  bT^Oiarriott  using 


a  T  - uy  ivieirriott  usinc 

a  cose  coupling  analysis  for  the  excitation  of  CO  Co'  enrt  u''rt  ■  . 

CO..  »'.t  u.  .  _  CO^  and  H  O  in  collision  with 


rr>  u  -h  tV  !  ,  '  noVvi  in  coiiisio 

CO^,  H^O,  and  He  are  reviewed  and  presented  in  a  form  convenient  tor 

use  in  modeling  the  chemistry  of  rarefied  flow  environments.  Although  no 
direct  cross  section  data  is  available  for  any  of  the  excitation  processes 
considered,  an  indirect  evaluation  of  the  theoretical  values  is  obtained  by 
comparison  to  available  relaxation  time  measurements.  .  I 


DD  Ii7a 

1  NOV  63 


Unclassified 

S«curity  ClaBsification 


(OVER) 


Unclassified 

Security  Cl aasiH cation 


ROLE  I  WT  I  HOLE  I  WT,  f  |  ROLE  I  WT 


Energy , transfer  cross  sections 
Vibrational  relaxation 


Unclassified 


Security  ClaRsij^Sdon 


A  REVIEW  OF  VIBRATIONAL  CROSS  SECTIONS 
INVOLVING  H2O,  CO2,  H2,  CO  and  He 


Edward  R,  ^'sher 


Wayne  State  University 
Research  Institute  for  Engineering  Sciences 
Departmert  of  Chemical  Engineering 
Detroit,  Michigan  48202 

Contract  No.  P  19628-72-C-0007 


SPECIAL  REPORT 
August  1972 

Contract  Monitor;  Alva  T.  Stair 
Optical  Physics  Laboratory 


Approved  for  public  release;  distribution  unlimited 


Sponsored  by . 

Defense  Advanced  Research  Projects  Agency 
ARPA  Order  No.  1856 
Monitored  by 

Air  Force  Cambridge  Research  Laboratories 
Air  Force  Systems  Command 
United  States  Air  Force 
Bedford,  Massachusetts  01730 


Introduction 


I . 

Under  conditions  in  which  the  molecular  velocity  distribution 
is  non-Maxwellian,  the  characterization  of  chemical  processes 
cannot  be  made  using  rate  constants  but  must  be  calculated  from 
the  fundamental  inelastic  cross  sections.  Since  these  non¬ 
equilibrium  environments  are  not  described  by  a  "temperature", 
the  excitation  produced  in  each  collision  must  be  separated, 
evaluated  and  summed  in  order  to  evaluate  the  total  excitation  rates. 
In  particular,  for  the  purpose  of  estimating  the  infrared 
radiative  emission  from  a  collisionally  excited  gas  when  the 
relative  collision  energy  is  high  (of  the  order  of  .25  to  several 
eV)  and  the  effective  rotational  temperature  of  the  colliding 
molecules  is  of  the  order  of  a  few  hundred  degrees,  a  cross 
section  description  is  required  to  calculate  the  initial  vibra¬ 
tional  excitation.  At  preser.-*-  no  experimental  procedure  has 
been  publicized  which  will  measure  vibrational  cross  sections  for 
neutral  species  of  interest  in  the  upper  atmosphere  polluted 
with  rocket  exhaust  effluents  (CO2/  ^2©,  H2,  CO,  0,  N2,  NO,  etc.). 
Without  this  experimental  input,  theoretical  methods  must  be 
relied  upon  to  furnish  excitation  cross  sections  for  relevant 
vibrational  modes  of  these  atmospheric  species. 

Over  the  past  10  years,  Marriott has  applied  a  close 
coupling  analysis  to  calculate  the  vibration-translation  (VT) 
cross  sections  for  many  of  the  internal  modes  of  CO2,  H2O,  CO 
and  H2  of  interest  in  high  altitude  plume  studies.  This  report 
summarizes  and  discusses  this  past  work  and  presents  some  new 
theoretical  cross  section  data.  In  order  to  provide  some 
experimental  verification  for  the  cross  section  calculations,  since 
no  direct  cross  section  measurements  are  available,  the  cross 
sections  have  been  integrated  over  a  Boltzmann  distribution  of 
relative  energies  and  the  resulting  rate  coefficients  compared 
to  available  rate  or  relaxation  time  measurements.  This  procedure 
is  not  entirely  satisfactory  since  the  current  close  coupling 
code  accounts  only  for  strong  VT  interactions  while  the  thermal! zed 
energy  transfer  processes  associated  with  species  like  H2O  and  H2 
are  strongly  effected  by  vibration-rotation  (VR)  and,  in  some 
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^  cases,  vibration-vibration  (W)  processes.  Thus,  the  calculated 

j  rate  coefficients  are  expected  to  be  no  more  than  a  lower 

I  limit  to  the  rate  coefficients  with  increased  agreement  expected 

j  at  higher  temperatures  where  VT  processes  become  more  important 

j  due  to  the  strong  temperature  dependence.  This  trend  is  found 

I  in  all  cases  considering  in  this  report  except  for  the  CO  +  h 

energy  transfer  system.  For  this  system,  it  would  ippear  that^ 
the  repulsive  intermolecular  potential  assumed  in  the  close 
coupling  calculation  is  significantly  too  steep.  A  further 
discussion  of  this  case  is  given  in  the  results  section.  The 
trend  toward  agreement  between  the  calculated  and  measured  rate 
coefficients  at  higher  temperatures  for  the  other  transfer 
systems  considered  in  this  report  suggests  that  the  calculated 
VT  cross  sections  are  reasonably  reliable  for  predicting  the 
vibrational  excitation  of  exhaust  species  by  atmospheric 
constituents  under  conditions  of  high  relative  collision  energies 
and  low  effective  rotational  temperatures. 

In  the  next  section,  a  brief  account  is  given  of  the  close 
coupling  analysis  and  of  the  analysis  of  the  theoretical  cross 
sections  in  terms  of  rate  coefficients  and  relaxation  times. 
Section  III  presents  the  results  of  the  comparison  of  calculated 
rate  coefficients  with  available  relaxation  time  measurements 
and  a  discussion  of  each  energy  transfer  system  for  which 
calculated  cross  section  information  is  available.  Lastly  in 
Section  IV,  a  summary  discussion  is  given  of  the  range  of 
reliability  of  the  calculated  cross  sectionr  based  upon  the 
calculated  rate  coefficients  and  evaluation  of  the  approximations 
used  in  the  close  coupling  analysis. 

A  complete  tabulation  of  the  excitation  and  de-excitation 
rate  coefficients  generated  from  each  cross  section  presented  in 
this  report  for  the  temperature  range  300  to  5000  K  is  available 
directly  from  the  author. 
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I i •  Theoretical  Analysis 

Since  complete  details  of  the  close  coupling  calculation  are 
available  elsewhere ,  only  a  brief  review  of  the  method  used 
v/ill  be  given  here.  The  excitation  collision  is  envisioned  as 
occurring  between  a  structureless  incident  molecule  and  a  molecule 
possessing  internal  structure  characterized  by  its  normal 
internal  modes.  The  coupling  between  rotational  and  vibrational 
states  is  neglected  and  the  scattering  potential  is  assumed 
to  be  spherically  symmetric  and  to  be  separable  in  terms  of  the 
vibrational  and  rotational  coordinates,  i.e.  the  breathing 
sphere  model  of  Herzfeld,  Schv/artz  and  Slawsky^^^.  For  most  of 
the  cases  treated  in  this  report,  the  intermolecular  potential  is 
assumed  to  be  given  by  the  Lennard-Jones  or  Stockmeyer  potentials . 
The  intramolecular  potentials  are  assumed  to  be  given  by  an 
exponentially  repulsive  form  in  which  the  exponential  coefficient 
is  given  by  the  shape  of  the  intermolecular  potential  at  the 
classical  turning  point The  vibrational  wave  functions  are 
taken  to  be  harmonic  for  each  independent  mode  of  vibration. 

The  total  cross  sections  are  obtained  by  summing  over  all 
significant  partial  waves.  In  general,  for  molecular  collisions, 
several  hundred  partial  waves  must  be  taken  into  account.  However, 
it  has  been  found  that  the  inelastic  cross  sections  normally  are  slowly 
varying  functions  of  the  partial  wave  number  and  relatively  few 
partial  cross  sections  are  required  to  permit  graphical  integration 
to  be  employed.  The  details  of  the  close  coupling  equations  and 
their  method  of  solution  are  reproduced  elsewhere . 

The  vibrational  state  designation  of  the  polyatomic  species 
considered  in  this  report  are  given  in  Table  I  and  the  energies 

associated  with  the  various  vibrational  modes  are  given  in 
Table  II. 
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Table  I. 

Vibrational  State  Designation  for  Polyatomic  Species 


Designation 


0 

1 

2 

3 

4 


Quantum  Number 


0  0  0 

0  1  0 

0  2  0 

loo 

0  0  1 
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Table  II. 

Vibrational  State  Energies* 


0 

1 

2 

3 

4 

C02 

0 

0.083 

0.166 

0.164 

0.291 

0 

0.198 

0.395 

0.441 

0.465 

CO 

0 

0.266 

0.532 

energies 

in 

eV '  s  above 

the  0 

state. 
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The  direct  resultjs  of  these  close  coupling  calculations 
produce  cross  sections  which,  for  comparison  with  available 
experimental  data,  must  be  converted  to  Maxwellian  rate  coefficients. 

To  accomplish  this  we  note  that  the  rate  coefficients  for  excitation 
of  the  s  to  t  vibrational  transition  is  given  by^^^^ 

00 

J^st  “  jf  vQg^(v)exp(-mvV2kT)  v^dv, 

where  m  =  reduced  collision  mass 

V  =  relative  collision  speed 
g^  =  degeneracy  of  the  t^"  state 

and  Qg^(v)  =  cross  section  for  the  s  to  t  process. 

Detailed  balancing  on  this  process  yields  the  relationship 

(2)  k^gg^exp(-E^/kT)  *  k^^g^exp (-E^/kT) . 

Experimentally,  the  relaxation  time  is  a  more  frequently 

measured  quantity.  This  relaxation  time  can  be  related  to  the  excitati 
rate  coefficient  through  an  harmonic  oscillator  assumption,  yielding 

(3)  PTg^  =  [1  -  exp  [-(E^-Eg)/kT]]  ^ 

where  is  the  relaxation  time  in  atm-sec  and  R  is  the  gas 

constant. 

In  order  to  generate  rate  coefficients  and  relaxation  times 
from  these  equations,  the  cross  section  versus  energy  values 
must  be  fit  to  an  analytical  function  and  integrated  over  a 
Boltzmann  distribution  of  relative  energies.  The  functional 
form  which  appears  consistent  with  the  theoretical  results 
is 
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(4) 


=  C3  *[An(E/E^)P 


* 


C4 


where  E  —  relative  collision  energy 
E^  =  threshold  energy 
and  Qg^(E)  =  cress  section  in  cm^. 

Eitting  this  function  to  the  available  cross  section  data  points 
has  led  to  Table  III. 

For  convenience  in  understanding  the  magnitude  of  these 
inelastic  cross  sections,  the  probability  per  collision  per  a 

given  transitive  Pg^(v)  can  be  generated.  This  quantity  is  de¬ 
fined  by 


(5) 


“el 

i 


Where  is  the  elastic  scattering  cross  section  and  the 

summation  is  over  all  available  inelastic  processes  from  the 
initial  state  s. 

Application  of  equation  (5)  together  with  the  inelastic 
cross  section  data  shown  in  Table  III  and  elastic  scattering 
cross  sections  taken  from  reference  10  leads  to  the  excitation 
probabilities  shown  in  Figure  1.  Note  that  only  the  probability 
for  exciting  the  lowest  vibrational  level  is  shown  in  the  Figure. 

As  discussed  in  the  introduction,  there  are  no  cross  section 
measurements  with  which  to  compare  the  calculated  results. 
Therefore,  the  derived  rate  coefficients  can  only  be  compared  to 
available  experimental  data.  The  specific  reactions  for  which 
data  exist  are  presented  in  the  following  section. 
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Table  III 


•  •  2 

Coefficients  to  Cross  Section  (cm  )  vs.  Energy  (ev)  Function 


REACTION 

ET  (ev) 

C2 

C3 

C4 

N2 

+  AR  0(01) 

2.89E-01 

1.093E01 

l.llOE-23 

0.0 

CO 

+  HE  0(01) 

2.660E-01 

3.G56E00 

3.263E-20 

-2.147E00 

CO 

+  HE  0(02) 

5.310E-01 

3.659E00 

8.189E-23 

-4.8^7E00 

CO 

+  HE  0(12) 

2.660E-01 

3.779E00 

1.299E-19 

-8. 742E-01 

CO 

+  H20  0(01) 

2.660E-01 

1.306E01 

1.449E-18 

3.183E-00 

CO 

+  H20  0(02) 

5.310E-01 

1.175E01 

2.189E-20 

1.764E-00 

CO 

+  H20  0(12) 

2.660E-01 

6.830E00 

4.082E-22 

-2.357E-00 

H20 

+  H2  0(01) 

1.980E-01 

3.937E00 

3.755E-16 

1.380E-00 

H20 

+  H2  0(02) 

3.950E-01 

4.387E00 

1.630E-16 

1.734E-00 

H20 

+  H2  0(03) 

4.410E-01 

2.628E00 

2.205E-18 

-3.184E-01 

H20 

+  H2  0(04) 

4.650E-01 

3.488E00 

1.183E-18 

4.102E-01 

H20 

+  H2  0(12) 

1.980E-01 

2. 398E00 

1.939E-16 

7.175E-01 

H20 

+  H2  0(13) 

2.430E-01 

4.588E00 

5.263E-17 

1.374E-00 

H20 

+  H2  0(14) 

2.680E-01 

5.631E00 

2.619E-17 

1.999E-00 

H20 

+  H2  0(23) 

4.560E-02 

3.992E00 

5.143E-18 

3.068E-01 

H20 

+  H2  0(24) 

7.000E-02 

4.180E00 

3.589E-18 

5.583E-01 

H20 

+  H2  0(34) 

2.440E-02 

9.229E00 

5.393E-18 

2.230E-00 

H20 

+  H20  0(01) 

1.980E-01 

5.160E00 

2.754E-17 

1.287E-00 

H20 

+  H20  0(02) 

3.950E-01 

3.355EOO 

8.759E-19 

-3.844E-01 

H20 

+  H20  0(03) 

4.410E-01 

3.260E00 

5.112E-20 

1.400E-00 

H20 

+  H20  0(04) 

4.650E-01 

3.134E00 

2.065E-20 

l.lOOE-00 

H20 

+  H20  0(12) 

1.980E-01 

3.682E01 

1.676S-18 

-1.709E-00 

H20 

+  H20  0(13) 

2.430E-01 

3.977E00 

2.756E-17 

1.144E-00 

H20 

+  H20  0(14) 

r.680E-01 

3.728E00 

5.328E-18 

9.232E-01 

H20 

+  H20  0(23) 

4.560E-02 

6. 307E00 

1.993E-16 

2.017E-00 

H20 

+  H20  0(24) 

7.000E-02 

5.153E00 

1.067E-16 

1.835E-00 

H20 

+  H20  0(34) 

2.440E-02 

1.157E00 

1.419E-16 

2.017E-01 

H20 

+  C02  0(01) 

1.980E-01 

1.098E01 

1.708E-19 

1.569E-00 

H20 

+  C02  0(02) 

3.950E-01 

8.663E00 

9.932E-21 

6.524E-02 

H20 

+  C02  0(03) 

4.410E-01 

1.481E01 

1.103E-19 

4.090E-00 

H20 

+  C02  0(04) 

4.650E-01 

5.116E00 

1.689E-23 

3.680E-00 

H20 

+  C02  0(12) 

1.980E-01 

4.923E00 

3.386E-19 

3.393E-01 

H20 

+  C02  0(13) 

2.430E-01 

2.119E01 

1.793E-16 

8.457E-00 

H20 

+  ':02  0(14) 

2.680E-01 

1.451E01 

9.282E-18 

5.192E-00 

H20 

+  C02  0(23) 

4.560E-02 

8.291E00 

2.392E-17 

1.948E-00 

H20 

+  C02  0(24) 

7.000E-02 

7.518E00 

5.091E-17 

2.180E-00 

H20 

+  C02  0(34) 

2.440E-02 

4.117E00 

1.289E-16 

9.704E-01 

C02 

+  H2  0(01) 

8.280E-02 

3.748E00 

4.420E-16 

1.244E-00 

C02 

+  H2  0(02) 

1.660E-01 

2.133E00 

1.142E-16 

3.772E-01 

C02 

+  H2  0(03) 

1.640E-01 

4.145E00 

9.133E-17 

1.176E-00 

C02 

+  H2  0(04) 

2.910E-01 

5.473E00 

2.600E-17 

1.721E-00 
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Table  con ' t 


REACTION 


ET  (ev)  C2 


C3 


C4 


C02 

+ 

H2 

Q(12) 

8.280E-02 

C02 

+ 

H2 

Q(13) 

8.090E-02 

Co2 

+ 

H2 

Q(14) 

2.090E-01 

C02 

+ 

H2 

Q(32) 

2.000E-03 

C02 

+ 

H2 

Q(24) 

l.;'60E-01 

C02 

+ 

H2 

Q(34) 

1.280E-01 

C02 

+ 

H20 

Q(01) 

S.280E-02 

C02 

+ 

H20 

Q(02) 

1.660E-01 

C02 

+ 

H20 

Q(03) 

1.640E-01 

C02 

+ 

H20 

Q(04) 

2.910E-01 

C02 

+ 

H20 

0(12) 

8.280E-02 

C02 

+ 

H20 

0(13) 

8.090E-02 

C02 

+ 

H20 

0(14) 

2.090E-01 

C02 

+ 

H20 

0(32) 

2.000E-03 

C02 

+ 

H20 

0(24) 

1.260E-01 

C02 

+ 

H20 

0(34) 

1.280E-01 

C02 

+ 

C02 

0(01) 

8.28CE-02 

C02 

+ 

C02 

0(02) 

1.660E-01 

C02 

+ 

C02 

0(03) 

1.640E-01 

C02 

+ 

CO  2 

0(04) 

2.910E-01 

C02 

+ 

C02 

0(12) 

8.280E-02 

C02 

+ 

C02 

0(13) 

8.090E-02 

C02 

+ 

C02 

0(14) 

2.090E-01 

C02 

+ 

C02 

0(32) 

2.000E-03 

C02 

+ 

C02 

0(24) 

1.260E-01 

C02 

+ 

C02 

0(34) 

1.280E-01 

2.267E00 

1.395E00 

3.476E00 

5.107E01 

7.104E00 

4. 391E00 

4.693EC'0 

5. 2 3 "200 

4.624E00 

4.696E00 

6.947E00 

0.495E00 

1.164E01 

1.030E01 

1.488E01 

4.601E00 

1.955E01 

1.144E01 

1.630E01 

1.926E01 

2.804E00 

1.064E01 

2.067E01 

1.048E00 

1.968E01 

8.858E00 


8.757E-16 

2.359E-17 

5.288E-18 

8.254E-38 

9.580E-18 

5.831E-18 

1.256E-16 

2.408E-17 

1.387E-18 

1.870E-20 

5.245E-17 

8.598E-18 

7.792E-16 

6.608E-21 

1.483E-15 

9.341E-19 

4.993E-19 

3.059E-20 

4.223E-20 

1.209E-17 

4.183E-18 

2.919E-19 

1.880E-15 

2.062E-17 

7.164E-16 

6.543E-20 


9.479E-01 
4. 906E-02 
4.379E-01 
7.115E-00 
1.683E-00 
3.820E-01 
1.173E-00 
1.217E-00 
2.679E-01 
-1.638E-00 
1.795E-00 
1. 726E-00 
5.282E-00 
1.440E-00 
6.078E-00 
-4.290E-01 
4.957E-00 
1.324E-00 
3.693E-00 
7.240E-00 
0.0 

2.411E-00 

8.969E-00 

0.0 

7.333E-00 
2. 766E-01 
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PROBABILITY  /  COLLISION 


COLLISION  ENERGY 


Figure  1:  vibrational  excitation  probabilities  per  collision 

as  a  function  of  collision  energy. 
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Ill .  Results 


1 


CO(v)  + 

« 

The  relaxation  time  measurements  of  Millikan  SsL;rve  as 
a  comparison  of  the  cross  section  data  shown  in  Table  III. 

The  temperature  dependence  of  the  calculated  relaxation  time 
as  derived  from  the  cross  sections  is  in  good  agreement  with  the 
experimental  value  f  r  all  temperatures  above  about  300  K  and 
follows^a  T  dependence  in  accordance  with  the  Landau-Teller 
theory  .  The  magnitude  of  the  relaxation  times  is  about  a 
factor  of  3  lower  than  experiment  indicating  that  the  cross 
sections  may  be  large  by  about  the  same  amount. 


This  rather  good  agreement  with  experimental  values  was 
obtained  after  the  use  of  experimentally  measured  interaction  - 
potential  data  '  as  opposed  to  the  Lenna.rd-Jones  or  Stock- 
meyer  form  used  throughout  most  of  the  other  calculations 
presented  here.  The  Lennard-Jones  results  gave  relaxation 
times  about  50  times  faster  than  experiment,  suggesting  that  the 
steeply  repulsive  r'^^  behavior  at  short  range  is  too  repulsive. 
The  improvement  by  using  the  more  weakly  repulsive  beam  data 
reflects  the  observed  fact  that  the  repulsive  potential  for  atoms 
interacting  with  molecules  is  usually  considerably  weaker  than 
r  (more  like  r~'  or  r”®)  . 


C0(v)  + 

Figure  2  shows  the  comparison  between  the  calculated 
relaxation  times  from  the  cross  section  data  and  available 
experimental  data  The  most  recent  data  of  von  Rosenberg, 

Bray  and  Pratt  obtained  at  shock  tube  temperatures  appears 
to  be  the  most  reliable.  Clearly  a  large  discrepancy  exists 
between  the  calculated  and  measured  results  particularly  at  low 
temperatures.  This  difference  is  undoubtedly  due,  at  least 
in  part,  to  the  neglect  or  VR  and  W  processes  in  calculating 


11 


PT  (ATM-SEC) 


T(*K) 


'SOOO  2000  1400  l<XX)  800 


•06  .07  -06  -09  -10  -11 


Figure  2;  Comparison  between  the  calculated  relaxation  time 
and  available  measurements  for  the  C0(v~)  +  H-0 
energy  transfer  reaction. 

]2 


the  relaxation  times.  Inclusion  of  these  processes  would 
tend  to  preferentially  decrease  the  lower  temperature  relaxation 
times  and  provide  a  better  agreement  with  data.  It  is  important 
to  note  that  the  extrapolated  high  temperature  relaxation  values 
are  in  fairly  good  agreement  with  the  extrapolated  data. 

At  high  temperatures  (above  about  400°-5000  ..a, 

-  cLr::::rr= 

collision  energies  above  about  0.5  ev. 


HoO(v)  + 

There  is  no  experimental  data  on  the  relaxation  of  H.O  by  H 
so  only  inferences  can  be  made  on  the  reliability  of  the  cross  ^ 
section  calculations.  Certainly  at  low  temperatures  (below  about 
2000  K)  calculated  rate  coefficients  from  the  cross  section  data 
are  likely  to  be  low  (relaxation  times  slow)  due  to  the  neglect 
of  VR  processes.  .Since  both  H2  and  H2O  have  low  moments  of 
ine.rtia,^VR  processes  are  likely  to  dominate  the  thermal  relaxation 
process  .  At  higher  temperatures,  VT  processes  become 

relatively  more  important  and  the  cross  section  data  should  become 
more  reliable. 


Also  based  on  the  C02(v2)  +  H2O  (H2)  relaxation  times  to 
be  presented  later  in  this  section,  the  high  temperature  calculation 
of  n^Ofv^)  +  H2  is  expected  to  be  pretty  good.  The  argument  for 
this  expected  agreement  is  somewhat  obtuse  and  assumes  that 
the  major  contribution  to  agreement  with  experimental  results 
comes  from  the  form  of  the  interaction  potential.  The 
argument  is  based  on  the  observation  that  the  calculated  CO 
{V2)  +  H2O  high  temperature  relaxation  times  are  in  good  ^ 
agreement  with  data  while  the  CO2  (V2)  +  H2  relaxation  times 
are  about  an  order  of  magnitude  faster  than  data  support 
(Vida  infra) .  By  contrast,  the  relaxation  times  are  about  an 
order  of  magnitude  too  slow.  Therefore,  since  the  inter¬ 
action  of  H2  and  H2O  with  CO2  and  H2O  are  expected  to  be 
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SOTiewhat  similar,  the  noted  discrepancies  in  the  CO.  and  H  0 
relaxation  results  would  be  expected  to  cancel  and  thus  the 
H20(v2'  +  H2  results  should  be  moderately  good.  The  expected 
reliability  should  be  better  at  higher  temperatures  where  VT 
processes  are  expected  to  be  important  in  the  relaxation  process, 
although  the  degeneracy  in  the  V2  mode  of  CO2  and  the  close 
resonance  existing  bet-.^een  the  mode  and  the  v  mode  of 
CO2  somewhat  weaken  this  agreement. 

H.,0(v)  + 

Tne  relaxation  of  H20{v2)  another  energy  transfer 

process  in  which  the  low  temperature  relaxation  behavior  is 
likely  to  be  dominated  by  VR  and  W  processes.  The  calculated 

considerably  longer  than  experimental 
values  .  This  discrepancy  is  clearly  shown  in  Figure  3  . 

As  in  the  earlier  comparisons  of  this  section,  the  agreement 

becomes  considerably  better  at  high  temperatures  where  VT  processes 
become  more  ir.portant. 

Unfortunately  due  to  the  extremely  fast  relaxation  times 
found  for  H2O,  high  temperature  data  has  not  yet  been  obtained. 
Based  on  a  simple  extrapolation  of  the  low  temperature  data,  the 
predicted  relaxation  times  are  about  an  order  of  magnitude  too 
slow.  This  implies  that  the  cross  sections  are  about  an  order  of 
nagnitude  too  small.  Until  additional  experimental  cata  become 
available,  a  more  reliable  assessment  of  the  calculated  cross 
sections  cannot  be  made. 

H^O(v)  + 

There  is  no  experimental  data  av-.ilable  on  this  relaxation 
process  but  some  qualitative  arguments  can  be  made  based  on  the 
reverse  relaxation  process,  i.e.,  C02(V2)+  H2O,  for  which  data 
exists.  As  will  be  shortly  demonstrated,  the  calculated  CO  (v  ) 

+  H2O  high  temperature  relaxation  times  are  in  good  agreement  ^ 
with  experiment.  This  implies  that  the  interaction  potential 
data  used  in  the  calculations  is  pretty  good.  Since  the  same 
interaction  potential  is  used  in  the  H20(V2)+  CO2  relaxation 
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Figure  3 :  Comparison  between  the  calculated  relaxation  time 
and  available  measurements  tor  the  H20(v2)  +  H  0 
energy  transfer  reaction. 
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calculation,  the  calculated  high  temperature  relaxation  times, 
and  hence  cross  sections,  are  expected  to  be  fairly  good. 

Again,  direct  experimental  data  in  the  high  temperature  region 
would  clarify  this  comparison. 

CO^(v)  + 

(21  22) 

The  available  experimental  data  on  this  relaxation  process  ' 

(see  Figure  4 )  suggest  that  the  calculated  relaxation  times  are 

about  an  order  of  magnitude  too  fast.  This  result  is  not  easy  to 

explain  since  VR  processes  are  expected  to  be  important  at 

lower  temperatures  which  would  make  the  disrgreement  even  greater. 

The  major  factor  in  thic  disagreement  appears  to  come  from  using 
-12 

the  r  Lennard-Jones  repulsive  potential  in  this  calculation. 

If  a  weaker  repulsive  force  law  were  used,  the  relaxation  times 
would  increase  and  better  agreement  would  result.  Thus,  from  this 
limited  comparison  the  calculated  cross  sections  should  be 
reduced  by  about  an  order  of  magnitude. 

CO^(v) 

Figure  5  shows  the  available  experimental  data  on  the  CO  (v«)+  H„C 

( 2  3—  26)  £  £> 

reaction'  '  in  comparison  to  the  VT  relaxation  calculations 

of  Marriott.  Again  VR  processes  involving  H2O  are  likely  to 

reduce  the  low  temperature  relaxation  time  calculations  of 

Marriott  bringing  better  agreement  with  data  below  about  1000  K. 

The  high  temperature  calculated  relaxation  times  appear  to  agree 

(23) 

with  the  extrapolated  data  of  Shields  and  Burke  and  with 

(24) 

Lewis  and  Lee '  '  but  lie  below  that  of  Bulthius  and  Ponsen 

Thus,  the  high  temperature  calculations  appear  to  be  fairly 
good  for  this  system.  The  cross  section  calculations  are  there¬ 
fore  expected  to  be  fairly  reliable  for  this  energy  transfer 

(24) 

reaction.  It  is  worth  noting  that  Bulthuis  and  Ponsen  show 
a  comparison  of  their  data  with  rate  coefficients  calculated 
from  Marriott's  cross  sections.  Their  calculation  based  on 
Marriott's  cross  sections  assume  a  linear  interpolation  between 
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Figure  4;  Comparison  between  the  calcular.ad  relaxation  time 
and  available  mearurements  for  the  C02(v2)  +  ^2 
energy  transfer  reaction. 
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Figure  5:  Comparison  between  the  calculated  relaxation  time 
and  available  measurements  for  the  C02(v2)  +  HjO 
energy  transfer  reaction. 


energy-cross  section  points.  The  formulation  used  in  this 
report  to  correlate  the  calculated  cross  section  data  is  expected 
to  be  considerably  better  than  a  linear  interpolation,  even 
though  the  rate  coefficients  calculated  from  this  cross  section 
function  give  slightly  poorer  agreement  with  the  data  of 
Bulthuis  and  Ponsen. 

CO^(V)  +  I 

A  large  amount  of  experimental  relaxation  time  data  Is 
available  on  the  C02(V2)+  CO2  relaxation  process  and  is  reviewed  in 
article  by  Taylor  and  Bitterman .  Figure  6  summarizes  this 
data  and  shows  the  comparison  to  the  calculated  relaxation 
times  of  Marriott.  Although  the  calculated  temperature  dependence 
is  somewhat  stronger  than  given  by  the  data,  the  magnitude  of 
the  agreement  is  pretty  good.  In  the  high  temperature  region 
the  calculated  relaxation  times  appear  too  fast  by  about  a 
factor  of  3.  Thus,  the  cross  section  data  should  be  reduced  by 
about  a  factor  of  3  above  0.5  eV  for  consistency. 
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Figure  6:  Comparison  between  the  calculated  relaxation  time 
and  available  measurements  for  the  C0»(v-)+  CO, 
energy  transfer  reactions. 
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Discussion  and  Summary 

In  the  preceeding  section,  relaxation  times  predicted  from 
the  calculated  VT  cross  suctions  are  compared  to  available 
relaxation  data.  This  comparison  is  only  made  since  no  direct 
cross  sections  have  been  measured  for  any  of  the  processes  * 

presented  in  this  report.  Although  on-going  experiments  at 
Cornell  Aeronautical  Laboratories  may  provide  valuable  information 
on  tne  magnitude  of  some  , f  the  Excitation  cross  sections 
considered  in  this  report.  The  cross  section  values  presented 
in  this  report  have  all  been  predicted  using  a  close  coupling 
code  based  on  a  structureless  particle  incident  on  a  spherically 
averaged  molecule.  Except  for  the  CO  +  He  t/alculation,  the  inter- 
molecular  potential  information  has  beei  estimated  from  Lennard-Jones 
or  Stockmeyer  potential  information.  The  parameters  for  these 
semi-empirical  potentials  are,  in  general,  determined  from 
transport  coefficients .  Thus,  the  potential  parameters  so 
determined  are  sensitive  to  a  different  collision  energy  range  than 
IS  of  interest  in  the  cross  section  calculations. 

In  addition,  no  direct  information  is  input  to  the  code  on  the 
shape  of  the  repulsive  intramolecular  potential  for  the  polyatomic 
molecules.  The  classical  turning  point  argument  used  by  Marriott 
and  earlier  by  Herzfeld  and  Litovitz^^^,  appears  to  be  the  only 
method  currently  available  fcr  computing  the  effect  of  different 
collision  partners  on  the  form  of  the  intramolecular  potential. 

Current  studies  are  aimed  at  either  clarifying  this  procedure  or 
suggesting  a  new  methodology  for  relating  intermolecular 
potential  data  to  the  internal  coordinates  of  the  molecule. 

As  pointed  out  in  the  theory  section,  the  interaction  is 
modeled  by  a  spherically  averaged  potential  weighted  by  the  method 
of  Keesom'  'for  preferential  collision  directions.  For  appli¬ 
cations  where  relative  collision  energies  are  large  (greater  than 
a  few  eV)  while  rotational  motion  is  characterized  by  low  tempera¬ 
tures  .  i.e.  effectively  a  molecular  beam  configuration,  this  potential 
averaging  procedure  should  be  replaced  by  calculating  cross  sections 
for  different  collision  angles  then  averaging  the  resulting  angle 
dependent  cross  sections.  For  the  collision  systems  of  concern  in 
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this  report  and  for  collision  energies  above  a  few  ev's,  it  is  not 
anticipated  that  strong  angular  dependent  interactions ■ are  dominant, 
but  further  studies  into  this  point  are  undejnvay.  It  should  be 
noted,  however,  that  the  largest  quenching  cross  sections  are 
invariably  obtained  for  non-zero  impact  parameters.  Yet  it  is 
only  for  zero  impact  parameter  that  a  fixed  orientation  can 
normally  be  maintained  during  a  collision.  Therefore,  some  type 
of  potential  averaging  procedure  will  alv.ays  be  the  most  realistic 
approximation  to  the  scattering  potential. 

In  general,  the  agreement  obtained  between  high  temperature 
rate  data  and  the  calculated  relaxation  times  is  rather  good. 
Although  indirect,  this  agreement  tends  to  support  the  hypothesis 
that  the  calculated  VT  cross  sections  given  in  Table  III  are  the 
most  rtllalle  source  of  excitation  data  at  the  present  time  for 
I’.se  in  the  chemical  modeling  of  rarefied  plume  environments. 

Current  studies  are  addressing  some  of  the  approximations 
used  in  the  codes  and  in  obtaining  0  atom  scattering  data  for  use 
in  calculating  the  vibratirnaJ.  excitation  cross  sections  for  plume 
species  by  atmospheric  0  atoms. 

A  complete  tabulation  of  the  calculated  vibrational 
excitation  rate  coefficients  over  the  temperature  range  300  to 
5000  K  is  available  by  writing  directly  to  the  author. 

Although  these  calculated  rate  coefficients  provide  the  only 
source  of  information  on  excitation  processes  involving  other 
than  the  mode  of  CO2  and  H2O,  the  predicted  values  are  subject 
to  considerable  uncertainty  due  to  the  neglect  of  both  VR  and  W 
processes.  Further,  scaling  of  the  higher  order  excitation  rate 
coefficients  by  the  discrepancy  found  in  the  V2  relaxation  time 
comparisons  may  be  considerably  in  error  but  this  procedure  may 
at  least  lead  to  order  of  magnitude  estimates. 
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